Majority-vote model on triangular, honeycomb and Kagome lattices 
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On Archimedean lattices, the Ising model exhibits spontaneous ordering. Three examples of these 
lattices of the majority- vote model with noise are considered and studied through extensive Monte 
Carlo simulations. The order/disorder phase transition is observed in this system. The calculated 
values of the critical noise parameter are qc = 0.089(5), qc = 0.078(3), and qc — 0.114(2) for 
honeycomb, Kagome and triangular lattices, respectively. The critical exponents P/v, 'y/u and 1/v 
for this model are 0.15(5), 1.64(5), and 0.87(5); 0.14(3), 1.64(3), and 0.86(6); 0.12(4), 1.59(5), and 
1.08(6) for honeycomb, Kagome and triangular lattices, respectively. These results differs from the 
usual Ising model results and the majority-vote model on so-far studied regular lattices or complex 
networks. The effective dimensionalities of the system Doff ~ 1.96(5) (honeycomb), Dcff = 1.92(4) 
(Kagome), and Doff ~ 1.83(5) (triangular) for these networks are just compatible to the embedding 
dimension two. 
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I. INTRODUCTION 

The majority-vote model (MVM) [l| defined on two- 
dimensional regular lattices shows second-order phase 
transition with critical exponents /3, ^, v — which char- 
acterize the system in the vicinity of the phase transition 
— identical [l|-l3| with those of equilibrium Ising model 

Bi- 

On the other hand MVM on the complex networks 
exhibit different behavior [6|-[l2|. Campos et al. investi- 
gated MVM on undirected small- world network [y] . This 
network was constructed using the square lattice (SL) 
by the rewiring procedure. Campos et al. found that the 
critical exponents j/v and /S/v are different from those of 
the Ising model [5[ and depend on the rewiring probabil- 
ity. Luz and Lima studied MVM on directed small- world 
network [7[ constructed using the same process described 
by Sanchez et al. |13[. They also found that the critical 
exponents j/iy and P/iy are different from these of the 
Ising model on square lattice, but contrary to results of 
Campos et al. Q for MVM the exponents do not depend 
on the rewiring probability. Pereira et al. |8| studied 
MVM on undirected Erdos-Renyi's (ERU) classical ran- 
dom graphs [ij], and Lima et al. 0] also studied this 
model on directed Erdos-Renyi's (ERD) and their results 
obtained for critical exponents agree with the results of 
Pereira et al. [8|, within the error bars. Lima et al. [lO| 
also studied this model on random Voronoy-Delaunay 
latti ce 1 1511 with periodic boundary conditions. Lima 
also [ll[ studied the MVM on directed Albert-Barabasi 
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(ABD) network 11 dl and contrary to the Ising model on 
these networks [I4], the order/disorder phase transition 
was observed in this system. However, the calculated /3/iy 
and j/v exponents for MVM on ABD and ABU networks 
are different from those for the Ising model [5| and de- 
pend on the mean value of connectivity z of ABD and 
ABU network. Lima and Malarz ^ studied the MVM 
on (3, 4, 6, 4) and {3^, 6) Archimedean lattices (AL). They 
remark that the critical exponents 7/1^, P/v and l/u for 
MVM on (3, 4, 6, 4) AL are different from the Ising model 
[5| and differ from those for so-far studied regular two- 
dimensional lattices [ll, 01 but for (3'',6) AL, the critical 
exponents are much closer to those known analytically 
for SL Ising model. 

The results presented in Refs. [6l-[l2| show that the 
MVM on various complex topologies belongs to differ- 
ent universality classes. Moreover, contrary for MVM on 
regular lattices [l|, |2| , the obtained critical exponents are 
different from those of the equilibrium Ising model [5|. 
Very recently, Yang and Kim [I3 showed that also for 
d-dimensional hypercube lattices (3 < d < 6) critical ex- 
ponents for MVM differ from those for SL Ising model. 
The same situation occurs on hyperbolic lattices |2C| . 

In this paper we study the MVM on three AL, 
namely on triangular (3^), honeycomb (6'^), and Kagome 
(3, 6,3, 6) lattices. 

The AL are vertex transitive graphs that can be em- 
bedded in a plane such that every face is a regular poly- 
gon. The AL are labeled according to the sizes of faces 
incident to a given vertex. The face sizes are sorted, 
starting from the face for which the list is the smallest 
in lexicographical order. In this way, the triangular lat- 
tice gets the name (3,3,3,3,3,3), abbreviated to (3^), 
honeycomb lattice is called {6^) and Kagome lattice is 
(3,6,3,6). Critical properties of these lattices were inves- 
tigated in terms of site percolation [2l| and Ising model 



Our main goal is to check the hypothesis of Grinstein et 
al. [23 — i.e., that non-equihbrium stochastic spin sys- 
tems with up-down symmetry fall in the universality class 
of the equilibrium Ising model — for systems in-between 
ordinary, regular lattices (like SL [l|) and complex spin 
systems (like spins on ERU and ERD [ll^ or ABU and 
ABD 0,111). 

With extensive Monte Carlo simulation we show that 
MVM on 6^, (3,6,3,6) and 3*^ AL exhibits second-order 
phase transitions with effective dimensionality D^s ~ 
1.96, 1.92 and 1.83 and has critical exponents that do 
not fall into universality class of the equilibrium Ising 
model. 



II. MODEL AND SIMULATION 

We consider the MVM [l[ defined by a set of "vot- 
ers" or spin variables a taking the values -|-1 or —1, 
situated on every node of the 6^, (3,6,3,6) and 3^ AL 
with N = L"^ sites for 6^ and 3^, and N = 3L'^ sites for 
(3,6,3,6). The evolution is governed by single spin-flip 
like dynamics with a probability Wi of i-th spin to flip is 
given by 



1 

Wi = — 
2 



1 - (1 - 2q)a^ ■ sign ^ aj 



(1) 



and the sum runs over the number z = 3 (for (3,6,3,6) 
and (6'^) lattices) and z = 6 (for (3^) lattice) of nearest 
neighbors of i-th spin. The control parameter < g < 1 
plays the role of the temperature in equilibrium systems 
and measures the probability of aligning against the ma- 
jority of neighbors. It means, that a given spin i adopts 
the majority sign of its neighbors with probability (1 — q) 
and the minority sign with probability q [l|, l8l-.12j| . 

To study the critical behavior of the model we define 
the variable m = X]i=i ^il^ ■ Ii^ particular, we are inter- 
ested in the magnetization Af , susceptibility x ^^nd the 
reduced fourth-order cumulant [/ 



M{q) ^ (|m|). 



X(g) ^ iV ((m^) - (m)2) 



V{q) ^ 1 - 



3(m2)2' 



(2a) 



(2b) 



(2c) 



where (• • • ) stands for a thermodynamics average. The 
results are averaged over the A^run independent simula- 
tions. 

These quantities are functions of the noise parameter 
q and obey the finite-size scaling relations 



X = L^/Vx(^), 



-7- = L ' fu{x), 
aq 



2 

(3b) 

(3c) 



where v, /3, and 7 are the usual critical exponents, 
fm,x,u{x) a-re the finite size scaling functions with 



X ^ {q ~ q,)L'/'' 



(3d) 



being the scaling variable. Therefore, from the size de- 
pendence of M and x we obtained the exponents /3/i> and 
7/1^, respectively. The maximum value of susceptibility 
also scales as L^/"^. Moreover, the value of q* for which 
X has a maximum is expected to scale with the system 
size as 



qc 



hL-^'" with 5 « 1. 



(4) 



Therefore, the relations (Pc]) and (JH) may be used to get 
the exponent 1/v. We evaluate also the effective dimen- 
sionality. Doff, from the hyper-scaling hypothesis 



2/3/1/ + 7/1/ = D, 



off- 



(5) 



A/ = L-^/-/™(x), 



(3a) 



We performed Monte Carlo simulation on the 6'^, 
(3, 6, 3, 6) and 3^ AL with various systems of size N = 64, 
256, 1024, 4096, 16384, and 65536 for 6^ and (3,6,3,6) 
AL and A^ = 192, 768, 3072, 12288, 49152, and 196608 
for 3^. It takes 2 x 10^ Monte Carlo steps (MCS) to make 
the system reach the steady state, and then the time av- 
erages are estimated over the next 2 x 10^ MCS. One 
MCS is accomplished after all the N spins are investi- 
gated whether they flip or not. The results are averaged 
over A'lun (20 < A^iun < 50) independent simulation runs 
for each lattice and for given set of parameters (g, A^). 



III. RESULTS AND DISCUSSION 

In Fig. [T] we show the dependence of the magnetiza- 
tion M , Binder cumulant C/4, and the susceptibility x on 
the noise parameter g, obtained from simulations on (6"^), 
(3, 6, 3, 6) and (3^) AL with A^ ranging from A^ = 64 to 
196608 sites. The shape of M{q), U, and x curve, for 
a given value of A^, suggests the presence of the second- 
order phase transition in the system. The phase transi- 
tion occurs at the value of the critical noise parameter 
qc- The critical noise parameter qc is estimated as the 
point where the curves for different system sizes N in- 
tercept each other [2J]. Then, we obtain qc = 0.089(5) 
and Ul = 0.578(3); qc = 0.078(3) and C/4* = 0.613(4); 
qc = 0.114(5) and U^ = 0.601(5) for 6^, (3, 6, 3, 6) and 3^ 
AL; respectively. 

In Fig. [3] we plot the dependence of the magnetization 
M* ~ M{qc) vs. the linear system size L. The slopes of 
curves correspond to the exponent ratio P/u according 
to Eq. ([5a|) . The obtained exponents are (3/i/ = 0.15(5), 



1.0 

0.8 
0.6 
0.4 
0.2 



0.0 
0.000 





(6^) 


" L=8 


s. 


- L=16 


\\ 


- L=32 


^\\ 


- L^64 


l\ \^^-.^ 


- L=128 


n\ N, ^~~-~ 


L=256 


lo^\^^~~~-- 


1 . 1 


V:^^^;^ 



0.035 



0.105 0.140 




0.0 
0.06 

(d) 

10^ 



0.08 



0.10 



0.12 



10' 



?^ 



10^ 



'i 



1 (6') 


L=8 

L=I6 


_ 


L=32 




L=64 


! /^^ 


L=12S 




^^^J,=256 


-rf^-^^ 


r 


1,1, 



1.0 

0.8 

0.6 

0.4 

0.2 

0.0 
0.054 

(b) 

0.7 

0.6 

0.4 

0.3 

0.1 

0.0 
0.059 

(e) 

lo' 



10' 



^ ....^^ 


(3,6,3,6) 


- — '^"^ ^^ 


^^ 


L=16 V 


~x^^^^ 


L=32 \ 


v^V^~^~~~""~--..> 


L=64 


l\\ ^^ ^"""^ 


L=I28 


\\^^^^"^--_ 


L=256 
1 


>=^^P= 



0.068 0.082 



0.096 



- 




^.^^ (3,6,3,6) 


: 


^ 


— 


L=8 


l\\\^^ 


— 


L=16 


l\ \ ^v ^^ 


— 


L=32 


1 \ \ \v 


— 


L=64 


1 \ \ ^\ 


— 


L=128 


\ V ^,^ ^~"^ 


■ 


L=256 
1 


VV^!!I:::--j_ 



0.071 



0.084 



0.097 




0.140 



0.13 0.14 



,085 



(g) 



0.090 0.095 



0.100 



10^ 



(h) 





(3,6,3,6) 


L=8 


' 




L=16 


- 




L=32 

L=64 


1 




L=128 




^^ 


L=256 


1 


^ 




— , r r 


1 







(3^) 


— 


L=8 
L=16 








— 


L=32 




~ 




— 


L=64 




^ 




=.^_ = 


L=128 
L=256 










. 



















1,1,1 



072 



0.078 0.084 



0.090 



0.11 
(i) 



0.115 0.12 



0.125 



FIG. 1: (Color on-line). The magnetization M, Binder cumulant U4,, and susceptibility x as a function of the noise parameter 
q, for 1/ = 8, 16, 32, 64, 128, and 256 lattice sizes, and with N = L sites for 6 (first column), 3 (second column) and N = 3L 
sites for (3,6,3,6) AL (third column). 





FIG. 2: Plot In \qc{L) — qd versus the linear system size L for 
6^ (circles), (3,6,3,6) (squares), 3® (diamonds). 



FIG. 3: Plot the dependence of the magnetization M* 
M{qc) vs. the linear system size L. 



0.14(3), and 0.12(4), respectively for 6^, (3,6,3,6) and 
3*^ AL. 

The exponents ratio 7/1/ at q^ are obtained from the 
slopes of the straight lines with -f/i/ — 1.64(5) for 6^, 
7/1/ = 1.64(3) for (3,6,3,6), and 7/1/ = 1.59(5) for 3^, 



as presented in Fig. 21 The exponents ratio •'y/iy at 
Qxn^aAN) are 7/1/ = 1.66(8) for 6^, 7/1/ = 1.62(5) for 
(3, 6, 3, 6), and 7/1^ = 1.64(1) for 3^, as presented in Fig. 

m 

To obtain the critical exponent l/iy, we used the seal- 




FIG. 4: Susceptibility at qc versus L for (6^), (3,6,3,6) and 
{3'') AL. 




FIG. 5: Susceptibility at g^^ax {N) versus L for 6^, (3, 6, 3, 6) 
and 3^ AL. 



ing relation Q . The calculated values of the exponents 
l/i/ arc l/f = 0.87(5) for 6^ (circles), l/v = 0.86(6) 
for (3,6,3,6) (squares), and l/i^ = 1.08(6) for 3^ (dia- 
monds) (see Fig. [2]). Eq. ([5]) yields effective dimension- 
ality of systems D^tf = 1.96(5) for 6^, D^s = 1.92(4) for 
(3, 6, 3, 6), and D^s = 1.83(5) for 3^ The MVM on those 
three AL has the effective dimensionality close to two 
contrary to ER classical random graphs (0.99 < Dcff < 
1.02) [|| or directed AB networks (0.998 < D^s < 1.018) 
[12 | with roughly the same nodes connectivity {z = 3) as 
for (6^) and (3, 6, 3, 6), and (z = 6) 3^ AL. 

The results of simulations are collected in Tab. H] 



IV. CONCLUSION 

We presented a very simple non-equilibrium MVM on 
6^, (3,6,3,6) and 3^ AL. On these lattices, the MVM 
shows a second-order phase transition. Our Monte Carlo 



simulations demonstrate that the effective dimensionality 
-Doff is close to two, i.e. that hyper-scaling may be valid. 
Finally, we remark that the critical exponents "f/iy, P /v 
and l/^/ for MVM on regular 6^, (3,6,3,6) and 3^ AL 
are similar to the MVM model on regular (3, 4, 6, 4) and 
(3**, 6) [13I and are different from the Ising model Q and 

TABLE I: Critical parameter, exponents and effective dimen- 
sion for MVM model on 6^, (3, 6, 3, 6) and 3'^. For complete- 
ness we cite data for SL Ising model as well. 





6^ 


(3,6,3,6) 


3« 


SL Ising 


qc 


0.089(5) 


0.078(2) 


0.114(5) 




P/u 


0.15(5) 


0.14(3) 


0.12(4) 


0.125 


7/^" 


1.64(5) 


1.64(3) 


1.59(5) 


1.75 


7/^' 


1.66(8) 


1.62(5) 


1.64(1) 


1.75 


1/u 


0.87(5) 


0.86(6) 


1.08(6) 


1 


Oeff^ 


1.96(5) 


1.92(4) 


1.83(5) 


2 



"obtained using x{N) at q = Qc 
obtained using x{^) a.t q = q* 
'^obtained using ratio 'y/u given by dependence x(^) ^t q - 



■■ qc 



differ from those for so-far studied regular lattices [l|, 01 
and for the directed and undirected ER random graphs 
IMlSi and for the directed and undirected AB networks [111 . 
Il2| . However, in the latter cases @, Sllll, Il2[ the scaling 
relations ([3]) must involve the number of sites iV instead 
of linear system size L as these networks in natural way 
do not posses such characteristic which allow for N oc 
L'^ {d G Z) dependence H^. For (6^), (3,6,3,6) and 
(3^) AL some critical exponents are much closer to those 
known analytically for square lattice Ising model, i.e. /3 = 
1/8 = 0.125, 7 = 7/4= 1.75 and ly = I, but except for 
v they differ for more than three numerically estimated 
uncertainties. 
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